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ABSTRACT 

Planetary systems consisting of one star and n planets with equal planet 
masses /i and scaled orbital separation are referred as EMS systems. They rep- 
resent an ideal model for planetary systems during the post-oligarchic evolution. 
Through the calculation of Lyapunov exponents, we study the boundary between 
chaotic and regular regions of EMS systems. We find that for n > 3, there does 
not exist a transition region in the initial separation space, whereas for = 2, 
a clear borderline occurs with relative separation ~ /i^/^ due to overlap of res- 
onances (Wisdom, 1980). This phenomenon is caused by the slow diffusion of 
velocity dispersion (~ t^/^, t is the time) in planetary systems with n > 3, which 
leads to chaotic motions at the time of roughly two orders of magnitude before 
the orbital crossing occurs. This result does not conflict with the existence of 
transition boundary in the full phase space of N-body systems. 

Subject headings: N-body Problem; Stability; Chaos; Diffusion. 



1. Introduction 

Planetary system is a paradigm of classical N-body problems. The pioneer work of 
Poincare([l|) indicated that a general N-body system with > 3 is not integrable. According 
to the work of Arnold(0), a typical near- integrable Hamiltonian system with more than 3 
degrees of freedom is topologically unstable, i.e, the action variables can diffuse away in the 
phase space during a sufficient long period. However, for a specific planetary system, when 
the instability is significant and how long the system can be stable are unsolved problems 
with great astronomical interests(0). 

One of the most interesting questions is whether there exists a transition region between 
regular and chaotic motions in a n-planet system. Gladman(0) found that for n = 2, the 
orbits of two planets with equal masses (/i) will be chaotic as long as their relative separation 
> 2/2^/"^. The scaling /i^/^ was first found by Wisdom([^) in a restricted three-body model 



2 



with a resonance overlap criterion. However, whether there exist similar relationships for 
systems with n > 3 remains unknown. 

In a previous paper((0), hereafter Paper I), we investigated the orbital crossing time 
and orbital diffusion of planetary systems with equal planet masses and scaled separation 
(we call them the EMS systems). We found that for EMS systems with n > 3, the velocity 
dispersion of the planets diffuses with time t as ~ t^^^. In this paper, we study the transition 
from regular to chaotic regions in the initial separation space of these systems. 



An EMS system consists of a star with unit mass (/iq = 1) and n planets with equal 
masses /i (;U <^ 1) moving around it. Assume all the planets are initially in circular and 
coplanar orbits with semi-major axes Oj, where i = l,...,n are sequenced from inner to 
outside. We put the innermost planet at ai = lAU, which is the average distance from the 
Sun to the Earth. The initial separation between the planetary orbits, denoted as fco, are 
equal when scaled by their mutual Hill's radii Rh, 



The initial mean anomaly Mj and longitude of perihelion zui of each planet orbit are chosen 
randomly. 

In an inertial coordinate with origin on the mass center of the n + 1 bodies ( the star 
and n planets), the Hamiltonian of the system is given as. 



where and pi = fiiTi {i = 0, ...n) are the position and momentum of the star and the plan- 
ets, respectively, = |rj — rj| is the distance between bodies i and j, G is the gravitational 
constant. 

The equations of motion corresponding to eq. ([2]) are integrated by a Runge-Kutta 
method with adaptive step-sizes. The orbits are simulated up to 10^ years, with an accu- 
mulated error of relative energy ~ 10^*^ for weak chaotic or regular orbits and ~ 10^^ for 
strong chaotic orbits in the case of n = 3 (Fig. la). The Lyapunov exponents are defined as 
the limit of x(t) when time t tends to infinity, with 
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Model and Method 
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Fig. 1. — Evolution of 3-plant EMS systems with mass — 10~ . (a) Variations of accumu- 
lated error of the relative energy with time. The three lines correspond to initial separation 
/co = 4, 5, 8, respectively, (b) Variations of x(^) for orbits with different initial separation Uq. 



where ^{t) and ,^(0) arc tangential vectors in the phase space at time t and 0, respectively, 
I • I is the norm of a vector. We solve the variational equations of eq. (2) along the solution 
to get Normalization of ^{t) at every fixed time interval is required to avoid overflow 

during the simulations. 



3. Numerical Results 

A notable phenomenon in an EMS system with n > 3 is for some orbits, S,{t) may reach 
a limit of medium value ( ~ 10"^ yr"^) at time t ~ 10^ yrs, then it grows to a large limit ( 
~ 10~^ yr~^ ) after a long period time ( ~ 10^ yrs), indicating that the orbits finally become 
quite chaotic after a long time in a transitional state (e.g., orbits with /cq = 6, 7 in Fig. lb). 
Moreover, we find that these transition orbits are in the border region between strong (with 
large Laypunov exponents) and weak (with medium ones) chaotic motions. 

The transition phenomenon of x{t) does not occur in EMS systems with n — 2. As shown 
in Fig. 2a, different integration time spans give almost identical border kc between chaotic 
and regular motions. For ji — 10~^, kc ~ 5.2, which roughly coincides with the location 
predicted by Wisdom and Gladman's criterion, i.e., kc — 2//^/'''/(2///3)^/^ 2.3//"^/^^ 5.5. 

The evolution of x(t) is linked with the velocity dispersion (a) in the planetary systems. 
As demonstrated in paper I, one of the major characteristics of planetary systems with 
n > 3 is that a diffuses slowly with ~ t^/'^ due to the interactions between planets. The 
growth of velocity dispersion (orbital eccentricities) of the planets enhances the planetary 
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Fig. 2. — (a) Variations of x{t) with ko at different time: 10^ yr (circles), 10^ yr (triangles) 
and 10^ yr (stars) in a 2-planet EMS system with mass // = 10~^. (b) Variations of x{t) with 
/co at t = 10^ yr (stars) in a 3-planet EMS system with mass /i = 10~®. The solid and dotted 
curves plot the functions — log(Tc) and 2 — log(Tc), respectively. 

perturbations, leads to the growth of x{t) after a transition stage, and finally orbital crossing 
may occur after the systems undergo a long period of strong chaotic motions. Define the 
orbital crossing time as the minimum duration that either orbit crossing ( > Uj for 
i < j) or close encounter {ry < Rh) occurs. In paper I, we derived an empirical formula of 
Tc for /co > 2.3, 10"^ < H< 10~^°. For planets in circular orbits, it gives, 

rji f 

log(^) = (-2 -0.27 log /x) + (18.7+ 1.1 log/x)log(-^). (4) 
yr 1.6 

Fig. 2b plots the x(lOVr) and quantities — log(Tc). A relationship log(x(10Vi')) ~ 2 — log(Tc) 
suggests that the time scale that an orbital becomes chaotic (i.e. Lyapunov time ) is roughly 
two orders of magnitude less than the orbital crossing time. The correlation indicates that 
the growth of x(i) is due to the slow orbital diffusion towards high velocity dispersion. Since 
growth of velocity dispersion is a robust procedure in a planetary system with n > 3, we 
believe that the transitional regular orbits will eventually become chaotic after sufficient long 
time, hence there is no clear boundary between the domains of chaotic and regular motions. 



4. Conclusion 

In this paper, we study the boundary between chaotic and regular motions in EMS 
systems with n planets. Our results indicate that the growth of velocity dispersion is robust 
as long as n > 3, which gradually turns the orbits initially in regular regions to chaotic 
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ones. Hence there does not exist a borderline in the initial separation space, which does not 
conflict with the existence of boundaries in the full phase space (position and velocity space) 
of N-body systems. 
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